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Abstract 

Munoz, Seidel, and Muga [Phys. Rev. A 79 , 012108 (2009)], following 
an earlier proposal by Poliak and Miller [Phys. Rev. Lett. 53 , 115 (1984)] 
in the context of a theory of a collinear chemical reaction, showed that 
suitable moments of a two-flux correlation function could be manipulated 
to yield expressions for the mean quantum dwell time and mean square 
quantum dwell time for a structureless particle scattering from a time- 
independent potential energy field between two parallel lines in a two- 
dimensional spacetime. The present work proposes a generalization to 
a charged, nonrelativistic particle scattering from a transient, spatially 
confined electromagnetic vector potential in four-dimensional spacetime. 
The geometry of the spacetime domain is that of the slab between a pair 
of parallel planes, in particular those defined by constant values of the 
third ( 2 ) spatial coordinate. The mean N th power, N = 1, 2, 3, . . ., of the 
quantum dwell time in the slab is given by an expression involving an Af- 
flux-correlation function. All these means are shown to be nonnegative. 
The N — 1 formula reduces to an 5-matrix result published previously [G. 
E. Hahne, J. Phys. A 36 , 7149 (2003)]; an explicit formula for N = 2 , and 
of the variance of the dwell time in terms of the 5-matrix, is worked out. A 
formula representing an incommensurability principle between variances 
of the output-minus-input flux of a pair of dynamical variables (such as 
the particle’s time flux and others) is derived. 


1 Introduction 

In 1984, Poliak and Miller [1] proposed a formula for computing the quantum 
mean time for the system of a forward or reverse collinear chemical reaction 
AB+C <-> A+BC to dwell between two spatial surfaces exterior to, and includ- 
ing, the reaction zone. Their formalism depended on calculating the algebraic 
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sum of weighted averages of the two-flux correlation function between the nor- 
mal fluxes of particles at the same, or opposite, pairs of the spatial surfaces, 
and was based on a microcanonical ensemble for choosing the density matrix 
of the quantum mixture. In 2009, Munoz, Seidel, and Muga ([2], Eq. (14)) re- 
formulated the two-flux method of Poliak and Miller to allow the estimation of 
mean powers of the dwell time in potential scattering in one spatial dimension; 
they in effect presumed a density matrix derived from a pure quantum state. 
They showed that the two-flux formalism could be manipulated and reduced to 
a known single-flux method for computing the mean dwell time, and that the 
mean-square dwell time formula in the context of Poliak and Miller’s approach 
was mathematically equivalent to a previously known, but distinct, formula also 
involving a two-flux correlation function. 

Munoz, et al. also advanced a formula for the mean-cube quantum dwell 
time in terms of Poliak and Miller’s two-flux correlation functions, and gave or 
implied another formula for the mean ./V th power of the dwell time in terms 
of IV-flux correlation functions. For N = 3, in the free-particle case, the two 
formulas gave different results. 

In the present work, a general formula for the mean iV th power of the dwell 
time for a particle between a pair of z=constant planes containing the interaction 
zone, in terms of IV-flux correlation functions, is proposed. The physical system 
is that of a charged, structureless particle moving in four-dimensional spacetime, 
and scattering from a specified, transient, spatially confined, electromagnetic 
four- vector potential. The formulas for N = 1 and N = 2 reduce to those of [2], 
and the results for all higher powers reduce to a straightforward generalization 
of [2], rhs of Eq. (18). An incommensurability principle will be derived that 
establishes a lower bound to the product of variances of the output-minus-input 
fluxes of a pair of distinct dynamical variables of the particle. 

The ensuing material is organized as follows: In Section 2, we shall define 
general IV-flux correlation functions for the probability four-currents at an or- 
dered set of N points, N = 1,2,3,..., of a simple charged particle moving in 
a time-dependent electromagnetic field. The correlation function will prove to 
have zero four-divergence with respect to each set of coordinates and associ- 
ated vector indices (generalizing [2], Eq. (9)). We shall show further that a 
weighted average of such a correlation function of normally-directed currents 
on the direct product of N three-dimensional spacetime boundaries equals, by 
an TV-fold application of Gauss’s theorem, to a 4IV-dimensional volume integral 
that is a weighted sum of the current correlations in the direct product of the 
interior spacetime volumes. Section 3 applies the method to propose a general 
formula for the mean IV th power of the dwell time; the formula entails the com- 
putation of the IV-flux correlation function of the particle’s probability current, 
and can be expressed either as an integral over the product of N copies of a 
four-dimensional slab-like spacetime domain, or as integrals over the product of 
N copies of the three-dimensional boundary of the domain. All mean powers 
of the dwell time calculated according to the general formula are shown to be 
nonnegative, thereby satisfying an important physical criterion. We then show 
that the general formula reduces to results given in [2] for the mean first and 
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mean second powers of the dwell time. In Section 4, an explicit formula for the 
mean-square dwell time in terms of the S-matrix is derived, and it is shown 
the mean-square dwell time is greater than, or equal to, the square of the mean 
dwell time — hence, the variance of the dwell time is well defined. There are 
four appendixes: In Appendix A, the properties of the U operators, which are 
expressed in the familiar form exp[— i(tb — t a )H/H] for time-independent Hamil- 
tonians, are studied particularly in the context of spacewise propagation of wave 
functions, and where the Hamiltonians are time-dependent. In Appendix B, a 
formula giving a lower bound to the product of the variances of two delocal- 
ized observables in terms of the mean value of a generalized commutator of the 
observables is derived. In Appendix C, an expression closely resembling the 
Aharonov-Bohm operator for arrival time is shown to arise from calculating the 
mean flux of time across a single 2 =constant surface in spacetime. Appendix D 
gives arguments that are intended to provide insight into the presence of flux 
correlations in the computations involving delocalized (i.e., for measurements 
not confined to a single timelike or spacelike hyperplane) observables. 

2 Single-particle, multiple-point-current corre- 
lation functions 

In this section we shall define, and study some properties of ./V th -orcler ( N = 
1,2,3,...) correlations of the probability flow currents of a simple charged par- 
ticle moving in a given background electromagnetic field. The space of quantum 
states will be that of suitably well-behaved functions of the four real coordinates 
labeling points in spacetime: 

£ = (d = (e°, e, e, a = (*, *, v, *0 = (*, a a a = (*, a a 

Matrix elements of operators between quantum states in Dirac notation will 
entail integrals over a four-dimensional domain D 4 of spacetime: 

($m ') D * = f d% f d 4 &$(£ 1 r<£i|fi|&>*(6), ( 2 ) 

Jd 4 Jd 4 

where we omit |fi on the lhs and substitute 5 4 (£i — £ 2 ) for (?i|^l? 2 ) on the rhs 
if 12 is the unit operator on spacetime functions. It will normally be the case 
that $(£) and \E'(^) in the above are both solutions of the Schrodinger equation. 
Convergence problems can be dealt with by restricting the operators used, the 
spacetime domain involved, or the wave functions. 

The dynamics is given by the nonrelativistic Schrodinger equation, and the 
components of the four-vector current operator have the usual form in terms 
of the derivative operators and the electromagnetic four- vector potential ([3], 
Eq. (24.39), and [4], Eqs. (11.16) and (11.22)). In fact, we have the specified 
four-vector electromagnetic-field potential 

(^(0) = (V(0,A(0) (3) 
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and the Schrodinger equation for *k(£) of a particle with mass in and charge e 


“£*<« = L 


£[ 


3 = 1 




HO s 

--^(0 *(0 

z or 3 c J 


( 4 ) 


We define the spacetime matrix of the four-current operator J7" M (!E!) at position 
(S) = (T, R) as follows: 


(f2,r 2 | l 7' i (T,R)|t 1 .r 1 ) = 5(h-T)S(T- tj 

( <5 3 ( ri — R)<5 3 (R — r 2 ) if fi = 0; 

x ^ 3 (r 2 -R)[^(4-ofc) 

( -S-A j (T, R)](5 3 (R-r 1 ), if M = j = 1,2,3, 


( 5 ) 


where the under-arrows indicate the direction in which the derivative acts when 
the operator is sandwiched into a matrix element. Note that the matrices of 
the J^(E) operators are all Hermitean, at least on the space of differentiable, 
complex- valued test functions on spacetime. One can also verify that, if both 
<&(£) and ’!'(£) are solutions of (4), then the four-divergence of the associated 
(and in general, complex-valued) current distribution is everywhere zero (the 
summation convention on a /.t = 0, 1,2,3 is operative henceforth): 

E_($|jm (S )|^) = o. (6) 

Now let \h(£;k), where k ranges over all real three- vectors, be a complete, 
orthonormal set of solutions to (4), in that at every time t, 


/d 3 ^(t,r 2 ;k)^, ri ;kr = <5 3 (r 2 - iq); 


( 7 ) 


orthonormality will be defined in terms of the matrix fF°(T) in a wave vector 
basis, 


(4'(k 2 )|J r0 (T)|4'(k 1 )) 

= j d 3 R(V( k 2 )|J°(T,R)|^(k 1 )) 

= 5 3 (k 2 -k 1 ). (8) 

An analogous flux-related matrix, called fF 3 (i? 3 ) for ^-propagation, will be de- 
fined in (All). 

The vector k will be taken to be the wave vector of an input plane wave 
such that the state \h(t, r;k) contains purely outgoing scattered waves. We 
can now define a U^ 1 ' 1 operator that propagates solutions of (4) in time, where 
the superscript (1) is intended to indicate that not free-particle states, but 
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interacting-particle states are involved; the symbol will be used for the 
corresponding operator in the interaction- free case — see (A7ab). We take 

(t 2 ,r 2 |i7 (1) |fi,ri) = J d 3 fc^(t 2 ,r 2 ;k)^(t 1 ,r 1 ;k)*, (9) 

which by (7) reduces to the unit operator on a time=constant slice when t± = t 2 . 
Note that the operator U ^ yields a Hermitean matrix in spacetime, in that 

<^ 2 , r 2 |t/ (1) |^i, ri)* = (ti,ri|[/ (1) |f 2 ,r 2 ). (10) 

We shall also make use of a generalized density operator which has a complete 
matrix in spacetime. Let p(k 2 ;ki) be a Hermitean, positive semidefinite ma- 
trix with unit trace that represents the particular mixture of plane-wave states 
that existed at times before the particle interacts significantly with the electro- 
magnetic field. The spacetime matrix of the density operator p (L is taken to 
be 

(t 2 ,r 2 |p (1) |ti,ri) = J d 3 ki J d 3 fc 2 T(f 2 , r 2 ; k 2 )p(k 2 ; ki)^(ii,ri; ki)*. (11) 

The matrix of both operators U ^ and p ^ evidently satisfy (4) on the left 
spacetime arguments, and the complex-conjugate equation applied to the right 
spacetime arguments. 

We now define a quantum average of an ordered N - point current correlation 
function: 


= TripWj^iZ^uWj^iZ^uW... J» n (Z n )}. (12) 

In view of (9) and (11), (12) reduces to 

(./"'•••"v(H i; ...;H a -)) /; = J d 3 k ■!••• J d 3 k N+1 p( kjv+i.ki) 

x (T(k 1 )|^(Si)|4/(k 2 )) • • • (Ik (kjv) | J'^ N (Hjv) ( kjv+i) ) • (13) 

Since in (13) each is sandwiched between two solutions of (4), (6) implies 

that the four-divergence of the lhs of (12) or (13) at each of the N points is 
zero: 

-^(Hij.^.jSjv)). = 0, for each s = 1,2, .. ,,N. (14) 

We suppose that the particle is effectively free on or near the boundary of, 
and everywhere outside, the slab S' 4 in spacetime, defined as follows: 

S' 4 = { (T, R) | — 00 < T,R\R 2 < 00 , z a < R 3 < z b }. (15) 

The boundary dS 4 of the slab comprises the two three-dimensional planes de- 
fined by 

dS 4 = {(T,R)| - 00 < T,R\R 2 < 00 , R 3 = z a or R 3 = z b }. (16) 


5 



We define the weighted current matrix as follows. Let w s (^ s ), s = 1,2, . . . , TV, 
be a sequence of suitable real- valued weight functions. Then we take 

(UW^(E s )\Ui) = (l/2)K(6)(6l^(S s )|6+i) 

+ (^\J^(E s )\^ s+1 }w s (Ui)], (17) 

where the symmetrization of the operator J with the weight function w is done 
because for n = 1, 2, 3, the current operator includes differentiations so generally 
doesn’t commute with w. We can use the delta-functions to replace the argu- 
ments of w by S; because of the particular bidirectionality of the differentiation 
process the terms involving differentiation of w cancel out, so that we can factor 
the w out of the quantum matrix element: 

(tf(k a )|W"*«(3 8 )|tf(k 8+1 )> = (Ik (k g ) \J ,ia (S s ) |\k (k s+1 ))tu s (S s ) . (18) 

Because of (6), the four-divergence of such a matrix is 

f) 8m 

^(*(k s )|W^(5 s )|*(k s+1 )> = (vk(k s )|^(S s )|T(k s+1 )^. (19) 

The weighted mean iV-order current is now well-defined as follows: 

N 

Wf-^(S 1 ;...;S Ar ) = (j«^(Hi;...;S,)).n«, s (S s ). (20) 

S = 1 

We further define the mean total divergence (creation) Wp of the quantity 
nfLi w s in the slab as the TV-fold integral over the slab of the TV-fold diver- 
gence of the current density W : 

Wp = f d 4 ^... f d*S N dN Wf -^(S i; . ■ B„). (21) 

Js 4 JS 4 '" 0:1 N 

In view of (19), we have 



An TV-fold application of Gauss’s theorem to (21) yields the following surface 
integral for Wp\ 



The latter is a weighted average over TV-flux correlation functions, generalizing 
the two- flux weighted averages used in [1] and [2]. 
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We remark that a sum of product weight functions will also lead to a sum 
of entities as in (21) and (22). The particular choices of weight functions for 
TV = 2, that is, — (1/2) |T\ — T 2 | and — (l/2)(Ti — T 2 ) 2 , provide alternative means 
of deriving [2], Eqs. (17) and (18), respectively, when expressed in the form of 
[5], Eq. (D2); note, however, that the first of these does not fit the pattern of a 
sum of products of functions of Tf and T 2 . 


3 Mean powers of the dwell time 


The methods of Sec. 2 permit the construction of expressions that are candidates 
for the expectation values of the iV th powers of the dwell time — see also the 
discussion in Sec. 5. We take the weight functions w s to be 

rc s (5 s ) = =T„ for a = 1 , ... , TV, ( 24 ) 

and factor these into (23), while calling the result Ts 4 **n- We have from (22) 

T s ***n = f ^Sr--- f d 4 E N (j°^° N (E 1 -,..-E N )) (25) 

Js 4 Js 4 p 

In view of the simple structure of the operator J°(S) (see (5), (25)), 

[ d 4 S(T(k)|J°(S)jT(k')) S 4 = f d 4 £T(?;krv]/(£;k') 

JS 4 J s 4 

= <tt(k)|*(k')) S 4 

= «*(k')|*(k)} S 4)*. (26) 

We now have 

Ts 4 **n = J d 3 k\ ■ ■ ■ J d 3 fc/v+ip(kjv+i; ki) 

x (4'(ki)|4'(k2))s4 • • • (4'(k7v)| v I'(kAr+i))s 4 (27) 

The cases N = 1 and TV = 2 are given or derived in [2], Eqs. (1) and (18), 
respectively. For larger TV, (27) reduces to, in the notation of [2], an expectation 
value of the operator (T D ) N . 

We shall now show that the rhs of (27) is real and nonnegative for each 
TV > 1. The proofs for even and odd TV differ slightly. For TV = 2 M = 2, 4, 6, . . ., 
we can rearrange the rhs of (27), as given by (13): 


Ts 4 ** 2 M = / d 3 kM + 1 


dr hi ■■■ drk M / drk M + 2 ■ d 3 k 2M +i 


x ^(T(k M+1 )|'F(k M+ 2))s4 • • • (’F(k 2 M)|’I'(k 2 M+i))s 4 ) 

x p(k 2 M+ijki) ^('F(ki)|4'(k 2 ))54 • • • (T(kM)| v I'(kM+i))s 4 ^ • 

(28) 
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For each kjvf+i the inner integral on the rhs of (26) is the expection value of a 
state (within large parentheses) with respect to the positive semidefinite density 
matrix p\ the result Ts4** 2 m is therefore real and nonnegative. 

For N = 2 M — 1 = 1, 3, 5, . . using (3), we rearrange as follows: 



x (^(£M;kM+l)*(^(kM+l)|^(kM+ 2 ))s 4 ' ' ' 
x ( v k(k2M-l)|'I'(k2M))s 4 ^ 

X p(k2M;ki)^(^(ki)|^(k 2 )) S 4 ■■ ■ 

x (^'(kM-l)|^'(kM))s 4 ^'(^M; kw)j • 


(29) 


For each fixed £m> the inner integral above is an expectation value with respect 
to p\ hence Tg4** 2M _ 1 is real and nonnegative. 

Classically, a particle’s trajectory can be contrived to be tangential to the 
boundary dS 4 , whence its mean dwell time and all powers would be zero; for the 
present system, we shall avoid tangentially-moving wave function components 
(fc 3 = 0) due to mathematical singularities there. Nevertheless, one could choose 
a p such that the particle moves fast enough normal to the slab that it spends an 
arbitrarily small, but necessarily nonzero, quantum time in S' 4 ; zero is therefore 
the greatest possible lower bound for the quantum dwell time for a slab of any 
thickness. 


4 Variance of the dwell time 

In order for there to exist incommensurability relations for the dwell time of a 
particle in S 4 and other dynamical variables, it is necessary at least that there 
exist a nonnegative variance V(T)g 4 of the dwell time, which should be given in 
terms of the mean square and mean dwell times as defined in (25), as follows: 

V{T)s* = r S 4„ 2 - (T S 4 Ml ) 2 . (30) 

The standard deviation A (T) of the dwell time for the domain S 4 will then be 

A(T) 5 4 = [V(T) S *] 1/2 . (31) 

We shall define presently an expression for the variance of the dwell time in S' 4 
that has the desired properties. The question of incommensurability of mea- 
suring the dwell time jointly with other dynamical variables is discussed in 
Appendix B. 



First, note that the result is trivial for the domain Q 4 defined below in (Dl), 
since 


^Q 4 **iV 



X 


(■ tb-t a ) N , for N = 1.2,.... 



Ti=tb 

Ti=t 0 


TN—tb 


(32) 


The variance V (T)q 4 is therefore zero. This simple result carries over only partly 
into the predicted averages of powers of (25 — 2 a ) (“normal travel distance”) of 
the end points of the slab S 4 . In particular, if there is total reflection of the 
particle approaching either boundary, the normal travel distance of the particle 
inside the slab is zero, while if there is transmission but never a reflected signal, 
the normal travel distance is (zb — z a ), with a continuum of possible outcomes 
for partial transmission and reflection; a more complicated result obtains for 
the mean of the square of the normal travel distance [6] . 

We now define the variance of the dwell time so that it is nonnegative and 
satisfies (30): 


V(T) S 4 = 


X 



J d 3 k 3 

(T(k 3 )|S°^(~ 2 )|T(k 2 )) 


x y5(k 2 ;ki) 


x 




T 5 4 Ml <5 3 (k 3 — k 2 ) 


<5 3 (ki — k 3 )T54** 1 J . 

(33) 


The variance so defined is nonnegative, as for each fixed k 3 , the remaining 
expression is an expectation value of a state-vector-like entity with respect to 
the positive semidefinite p, and multiplying out the integrand yields four terms 
that reduce to (30). 

Next we shall give the formulas for Ts***i and Tg 4** 2 in terms of the S- 
matrix. The former was obtained in [7], Eqs. (89)-(92) (the closed-channel 
terms of [7], Eq. (93) play no role here); inasmuch as [7] used the wave vector 
comprising (k*, k x , k v ) instead of the (k x ,k v , k z ) used herein, we shall state the 
old result in terms of the new variables. It is convenient first to introduce a 
Hermitean operator-function A (zb,z a ), which acts on suitable complex-valued 
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functions in wave- vector space, that appears in both formulas: 


Rl=z b 


<tt(k5)|A(z 6 ,* 0 )|tt(k4)> = / dT 3 dR\dRl^(^)\T 3 J 3 {T^)\^(^)) 

J dS 4 


= / d 3 k fi 


! (k 5 — kg) 


i id 
2 (k 3 ) 2 + kldk 3 


+ H-kl)z^-S(kl)z. y (ke _ kt) 

S' t (k 5 - k 6 ) * * d , 8(H)z b - 9(-kf)z a 


2{k 3 ) 2 k 3 dk 3 


+ A 3 (k 5 — kg) 


7,3 

^6 


S(k 6 ;k 4 ) 


2fc||fc|| 


x (d(—kg) exp(— - 0(k%) exp(-2iklz a )^S(kl, fcg, -fcf; k 4 ) 
+ 5t(k 5 ;^,^,-fc|)(-^3 

x (d(-kl) exp(2 ik%z b ) - 0(kg) exp(2ifc|z a )^(5 3 (k 6 - k 4 ) 


(34) 


Note that the above is a differential operator that will act on arguments of the 
density matrix in the formulas to follow. 

We remark that the terms linear in S or S ^ in (34) represent interference 
between forward- and backward-traveling waves at one or the other boundary 
plane. For an electron with energy E = 10 _3 eV, the rationalized wave number 
(2 mE/ti 2 ) 1 / 2 « 2 x 10 8 m -1 , so that the thickness (z b — z a ) of a slab would have 
to be roughly 10 nm or less in order that the exponents in these interference 
terms be small enough that oscillations — due to unavoidable slight smearing 
out of the energy distribution — do not average these terms to zero. 

We remark also that in the free-particle case (unit S'-matrix), the operator 
A (zb,z a ) reduces to (see also (C3)) 


(H/(k 5 ) | A(^b, z 0 )|'F(k 4 ))| free = 5 3 (k 5 -k 4 )(^3^) 

+ S\kl - k\)8\kt - kl)S\kl + k\) 
x 2fc3 | fc3 | (exp(2iklz b ) - exp(2iklz a )) . (35) 


We can now state the formulas for the mean and mean square dwell time in 
S 4 : 

T s *„i = jj d 3 k i J d 3 fc 2 (^(k 1 )|A(^, 2a )|^(k 2 ))p(k 2 ;k 1 ), (36) 

and 

T s 4**2 = (f) 2 / d3k i J d 3 k 2 J d 3 /c 3 (^(k 3 )|A(0 b ^ a )|^(k 1 ))* 

x (vF(k 3 )|A(^ b ,^ 0 )|vI/(k 2 )) / 5(k 2 ; k 4 ). (37) 
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In (37), the A-operators both differentiate only the density matrix, so their order 
is immaterial. It is worthy of note that the formula (37), and others (see App. B) 
herein, can involve the calculation of mean values of a vector-like entity, given 
the positive semidefinite density matrix p, where the “vector” can be a linear 
partial differential operator in k-space; given suitable differentiability properties 
of p(k 2 ;ki), it is plausible, but, so far as I can determine, not yet subject to 
mathematical investigation, whether the positive semidefinite property persists 
when the vector is of this type or possibly a generalized function. In the case 
of (36) and (37) the nonnegativity of these expressions was demonstrated in 
the mathematically distinct forms of (28) and (29), where the “vectors” are, or 
seem to be, ordinary functions. 


Appendix A: U operators in spatial propagation 

In this Appendix we shall elaborate on the properties of plane-wave states, 
and on the free-particle U^ 0 ' 1 and interacting particle U ^ operators in spatial 
propagation of quantum states. 

First we consider the plane-wave states, which are solutions of the free- 
particle Schrodinger equation. 

ip(t, r;k) = (27r) -3 / 2 exp[ik • r — iKk • kt/(2m)]. (Al) 

These functions are orthonormal and complete on t=constant slices of spacetime: 


/d 3 r^,r;k im ,r;k 2 ) = <5 3 (k! - k 2 ), (A2a) 

/^#,n;kWn 2 ;k)« = <5 3 (r ! — r 2 ). (A2b) 


The same plane waves are also orthogonal and qualifiedly normalized, but not 
complete, on r 3 =constant slices of spacetime, provided that we use the 3- 
component of the current operator instead of the unit operator as the entity 
with respect to which the matrix elements are calculated. We have 


J dtdr'drMtyyy-, 


h , 

f d 

2 im 

^<9r 3 



k 2 ), 


= 7r~{ k l + k%)5(k{ - kl)8(k\ - kl) 

Am 

X ~ ( fc 2 ) 2 )) exp [ir 3 {kl - kl)} 

= a(k 3 )8 3 ( ki-k 2 ) 


(A3a) 


where in the last step no contribution arises from the zero of the argument of 
the third (5-function at fc 3 = —kf, and 

a(k 3 ) = k 3 /\k 3 \. (A4) 
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As stated, we must avoid plane-wave input contributions from k 3 ss 0 in prob- 
lems involving S' 4 . The plane wave states are not complete on the r 3 =constant 
surface, however, as only states of positive energy (= hk°), which can be con- 
strued as the negative of the momentum conjugate to time, are involved. Since 
for “on shell” states 


fc 3 = ±[(2 mk°/h) - (fc 1 ) 2 - (k 2 ) 2 } 1/2 , (A5) 

those states for which the argument of the square root is negative have real 
exponential r 3 -dependence, i.e., are closed channels. Closed-channel contribu- 
tions appear in the Green’s functions, but not in the If 1 ' 0 ' 1 operators — see [5], 
Eqs. (A17) and (A21). 

The propagation operator U^J 3 gives rise to the following evolution of 
free-particle states in the z direction: 

I d 4 6 J d% I dTdR l dR 2 (b\uW\b) 

x (^ 2 \J 3 (T : R\R 2 ,R 3 )\( 3 )^ 3 ) 

= (n + ^)(a)-(n_^)(6), (A3b) 

where ?/> is any solution of the free-particle Schrodinger equation having only 
real wave-vector (i.e., open-channel) components, and II + and II are projection 
operators on the “positive k 3 ” and “negative fc 3 ” parts of its operand —this 
result follows at once from (A3a). The minus sign for free-particle states with 
negative k 3 is desirable physically in that, for all terms in, say, Tg 4**1 or T54 ** 2 , 
the particle passes the boundary planes in reverse order as those with positive 
k 3 , and dwell times must be nonnegative. The calculation is more complicated 
in the interacting case as both reflections and transmissions can occur. 

Let us next consider the operator [A 1 -* in the temporal evolution ( Q 4 ) envi- 
ronment. With a plane wave (Al) as the input state, we infer from [7], following 
a translation of notation, that, given an input wave vector ki , the S'-matrix plays 
the following role in the interacting-state wave functions of (7), (8): 

T(f, r;ki) = ip(t, r;ki), for t < t a , (A6a) 

T(f,r;ki) = f d 3 k 2 i>{t 1 r; k 2 ) 5'(k 2 ; ki), for t > tb- (A6b) 
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We can now find explicit forms for U i 1 ) outside or on DQ 4 : 

(t 2 ,r 2 \U (1) \t 1 ,ri) = (t 2 ,r 2 |17 (0) |ti,ri) = J d 3 kip{t 2 ,r 2 \ k)^(^i , ri ; k)*, 
for t‘i < t a and t\ < t a , or for t 2 > t-b and t\ > tb, (A7ab) 

(t 2 ,r 2 |[/ (1) |ti,ri) = J d 3 k 2 J d 3 h ip(t 2 , r 2 ; k 2 )S(k 2 ; ki)^(ti, iq; ki)*, 

for t 2 > tb and t,± < t a , (A7c) 

(t 2 ,r 2 |f/ (1) |ti,ri) = J d 3 k 2 j d 3 k\ ^(t 2 , r 2 ; k 2 )S f (k 2 ; ki)^(ii, n; ki)*, 

for t 2 < t a and t\ > tb, (A7d) 

where we used the fact that S is unitary. 

Let ki and k 2 represent input and output wave vectors, respectively. The 
matrix S'(k 2 ;ki) comprises four blocks, according to the algebraic signs of its 
arguments and k 3 , as follows: 

k 3 < 0, k 3 > 0 : output at r 3 = z a , input at r 3 = z a , (reflection), (A8a) 

k 3 > 0, k 3 < 0 : output at r 3 = Zb, input at r 3 = Zb, (reflection), (A8b) 

k 3 > 0, fc 3 > 0 : output at r 3 = Zb, input at r 3 = z a , (transmission), (A8c) 

k 3 < 0, k 3 < 0 : output at r 3 = z a , input at r 3 = Zb, (transmission). (A8d) 

(The above correspond to the open-channel reflection and transmission matrices 
of [7], Eqs. (73) and (74), with changes in variables and notation.) 

We shall now specify the wave functions T of (7), (8) outside or on dS 4 , in 
three steps: We define 

[ +1, if r 3 < z a , 

C(r 3 ) = < undefined, if z a < r 3 < Zb, (A9a) 

(-1, if r 3 >z b . 

We also define 

* ±1 (k 2 ; ki) = d(±k 3 )S 3 ( k 2 - k,) + 6(Tk 3 2 )S(k 2 ; k 1 ). (A9b) 

where 9 is the unit step function. The desired result is 

T(f,r;ki) = j d 2 k 2 ip(t, r; k 2 )X c(r 3)(k 2 ; ki). (A9c) 

We note that (A9) entails the assumption that so long as r 3 is outside of 
S' 4 , but t ventures inside Q 4 , the properties of the S-matrix guarantee that 
the the time development of the wave functions in (A9) is still given correctly; 
conventionally, the S-matrix is developed in terms of its time dependence, as 
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an input-to-output map valid only in mapping the pre-interaction to the post- 
interaction wave function, as in (A6b). We are presuming that the same S- 
matrix also maps wave functions correctly when the input and output surfaces 
are merely exterior to S 4 , whatever be the choice of time variable. 

Let us define and calculate the k-space matrix of the net probability flux 
R 3 (R 3 ) flowing across planes R 3 =constant [compare (8)], for R 3 < z a or for 
R 3 > Zb- We find, respectively, that 

(vh(k 2 )|^ 3 (i? 3 )|vl/(k 1 )) 

= j dTdR 1 dR 2 ('f(k 2 )\J 3 {T,R 1 ,R 2 ,R 3 )\'k(k 1 )) (AlOa) 

= C(i? 3 )^(C( J R 3 )fc 2 3 )^ 3 (k2-k 1 ) 

- J d 3 k 3 St(k 2 ; k 3 )0(-C(i? 3 )fcf )S(k 3 ; kr)j ; (AlOb) 

the equality of the two forms of (AlOb) is equivalent to the unitarity of the 
S'-matrix. In fact, because of (6), once the rhs of (AlOb) is evaluated at an 
R 3 < z a , say, these values hold for any R 3 , including z a < R 3 < Zb- When the 
^-matrix is block-diagonal (no reflections), the rhs of (AlOb) reduces to (A3a). 
If the ^-matrix is block-off-diagonal (no transmissions), the rhs of (AlOb) is 
zero, corresponding to the fact that there is then zero net current for all A 3 ; 
hence, the combined spatial propagation operator U^J/ 3 zeroes wave functions 
in systems for which no transmissions occur, and is of little use in these cases. 
That the latter outcome is not associated with zero for the powers of the dwell 
time arises from the fact that in carrying out the integrals in (23) with the 
weight functions (24), it is not the net probability current but the net time 
current that is calculated — what results is a differential operator, as in (35) . 

We now exhibit how the U C) operators of (A7) propagate the wave functions 
of (A9) along z, from R 3 to r 3 : 

J dTdR 1 dR 2 (t,r\U m J 3 (T,R 1 ,R 2 ,R 3 )\'S>{k 1 )) 

= j d 3 k 2 'k(t,r-k 2 )(^(k 2 )\R 3 (R 3 )\^(k 1 )). (All) 

The last result is independent of R 3 . 

Appendix B: Lower bound for 
incommensurability 

We study the question of mean values and variances of system dynamical 
variables analogous to the dwell time, and of the compatibility or incommensu- 
rability of joint measurements of pairs of such dynamical variables, one of which 
can be the dwell time. “Joint measurements” does not mean simultaneous, as 
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we are working primarily with spatial rather than temporal propagation, and 
the measurements determine flux correlations that occur at different times and 
spatial positions. In this connection, note that, in the Heisenberg picture, dy- 
namical variables that commute at equal times generally do not commute when 
taken at unequal times. For example, the difference in r 1 positions of (D2)(i), 
will normally not commute with the similarly defined difference in r 2 positions. 
In time intervals such that the interaction-free realm is attained by the particle 
both at input and at output, one might use the S'-matrix to compute standard 
deviations and quantitative estimates of incompatibility for the measurement 
pairs; instead, we shall study measurements across z intervals and the slab S' 4 . 

Let K and L be Hermitean operators for which we want to compute mean 
and mean square differences of their values upon exit from, and entry to, the 
slab S 4 ; we shall also define an expectation value for a generalized commutator, 
and show that the product of the variances of K and L is greater than, or equal 
to, one-fourth of the square of the mean of the commutator. 

Since the operators do not generally commute with J 3 , we use a symmetrized 
version to maintain the Hermitean property. We shall not construct a general 
current density for I\ in the interior of S 4 , but work with boundary fluxes only. 
The average net output minus input flux of K for S 4 we call Kqs 4 ** 1 i which is 
taken to be 


KdS 4 **l = 


dTdR 1 dR 2 


las 4 


x Tr ((1/2) (KJ 3 (E) +J 3 (Z)K)pW 


R —Zb 
R 3 —z a 


= / d 3 k i / d 3 k 2 


dTdR 1 dR 2 


IdS 4 


x (^(k 1 )|(l/2)(^J 3 (S) + J 3 (S)^)|^(k 2 )) / 5(k 2 ;k 1 )|^. (Bl) 
Similarly, we find for mean square of the flux difference Kgg 4** 2 


K dS 4** 2 = [ d 3 k 1 [ d 3 k 2 [ d 3 k 3 [ dT^RldR^ [ dT 2 dR\dR 2 2 

J J J JdS i Jos 4 

x (^(k 3 )|(l/2)(^J 3 (S 1 )+ J 3 (S 1 )i^)|M/(k 1 ))*|^ 

x (T(k 3 )|(l/2)(/T l 7 3 (S 2 ) + J 3 (Z 2 )K)\*(k 2 ))\^Z Z z b a 

xp(k 2 ;ki). (B2) 

The variance V(K)gg4 and standard deviation A(K)gg4 are calculated as in 
(30), (31), and (33), using the symmetrized operator product of K and J 3 
instead of the time operator. 

Next we follow, in part, a procedure given by Stoler and Newman [8] to 
obtain a qualitative estimate of the incommensurability of measurements of K 
and L. Let A be a real number, and define 

M( A) = K + iXL = M(- A) f , (B3) 
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where the f means Hermitean conjugate in the spacetime sense. We take the 
variance of M( A) to be 


V(M(X)) as * = J d 3 h J d 3 k 2 j d 3 k 3 

x [ dT 1 dR 1 1 dRi(^(k 3 )\(l/2)(M(X)J 3 (E 1 ) 

l JdS i 

+ J 3 (S 1 )M(A))|^(k 1 ))|^ - S 3 (k 3 - k 1 )M(X) dS , 
dT 2 dR 1 2 dR 2 2 (^(k 3 \(l/2)(M(X)J 3 (E 2 ) 


**1 


L JdS 4 


+ J 3 (E 2 )M( A))|^(k 2 )|^‘ - <5 3 (k 3 - k 2 )M(X)ds 4 **i 
x p(k 2 ;ki). 


(B4) 


Since for every fixed k 3 in the integrand, the remaining expression in (B4) 
is, in effect, the expectation value of a function in k-space with respect to p, 
V(M(X))qs 4 is nonnegative. If we use (B3) in (B4), multiply out, and combine 
terms, we find that 

V(M(X)) dS * = V{K ) as 4 + A 2 V{L) as < + XC(K, L) dS 4, (B5) 


where C(K, L)gs 4 is an expectation value of an entity that is plausibly described 
as a generalized commutator of K and L, adapted to dS 4 : 


C(K,L ) dS 4 = - [ d 3 k i [ d 3 k 2 [ d 3 k 3 [ dT^RldRj [ dT 2 dR\dR ! 
4 J J J Jos 4 JdS 4 

x { [<v&(k 3 )| ( KJ 3 (~ 2 ) + J 3 (S 2 )K) |4/(k 2 ))p(k 2 ; ki) 
x <^(k 1 )|(LJ 3 (S 1 )+J 3 (S 1 )L)|^(k 3 )) 


r i ' 

i i-Ri=*i> 

K^L 



R.2=Zb 

R 3 o=z a ' 


(B6) 


Note that C(K, L)ggi is real. The quadratic expression in A of (B5) is nonneg- 
ative for all A, so that the discriminant of the expression must be nonpositive, 
whereupon 

V{K) dS * V(L) dS * > (1/4) (C(K, L) dS i) 2 \ (B7) 

(B7) is the desired inequality that sets a lower bound to the product of the 
variances of two output-minus-input fluxes on dS 4 . 

Something like a true uncertainty principle for arrival, or crossing, times 
of the particle at the R 3 = Zb plane might arise as follows: Choose K and L 
to be smoothly independent in a way that they both take the value zero at 
R 3 = z a and have a nonzero value at R 3 = Zb, with zero i? 3 -derivatives at both 
z a and Zb- Then the values of (Bl), (B2) and (B6) all have zero contributions 
from i? 3 2 = z a . The problem then becomes that of choosing A", say, to be T 
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at R 3 = Zb, and L to be something like the conjugate momentum px also at 
R 3 = Zb, and trying to show that (B6) is nonzero for any physical p , including 
pure states. I have not attempted this line of investigation. 


Appendix C: Aharonov-Bohm operator 

In 1961, Aharonov and Bohm proposed an operator for the mean arrival 
time for a free particle at a specified (one dimensional) spatial position, given 
the quantum state at some initial time: see [9], [10], and for recent work, [11]. In 
this appendix, we state that what is in effect that operator can also be defined 
on single z=constant slices of spacetime. Let r(z)p be the mean value of T=time 
on the slice R 3 = z of spacetime. Then we have, omitting intermediate steps, 


r(z)p = j d 3 k\ J d 3 k 2 J dTdR 1 dR 2 

x {ip{ki)\T J 3 {T , R 1 , R 2 , R 3 )\tjj(k 2 )} 


R*=z 


m 


+ 5(k\ - k\)5{k\ - kl)8{k\ + fcf) 


= - / d 3 k i / d 3 k 2 


p(k 2 ;k!) 

i d 


L2fe||fc|| 1*11 dkl 


■ 2k 3\ k 3\ e x P (‘liklz)] }p(k 2 ;ki). 


(Cl) 

~l3 
Ap 2 j 


(C2) 


In (Cl), the -i/>(k)’s are the free-particle wave functions of (Al). We remark that 
(35) is related to (C2) by 


r(z b )p - t(Zo)p 


j j d 3 h J d 3 fc 2 (T(k 1 )|Ap b , 2a )|T(k 2 ))| free p(k 2 ;k 1 ). 

(C3) 


Appendix D: The basics of classical-to-quantum 
correlation functions 

In this section we shall attempt a physical motivation for the appearance of 
correlation functions in formulas for certain expectation values. 

Conventional matrix elements of quantum operators specify values of the 
transition amplitude between two quantum states at a given instant of time. The 
mean dwell time for a particle in, say, the slab S' 4 is a more delocalized entity 
that in effect involves integrals over a spacetime domain, or of the boundary 
of that domain — two disjoint, parallel three-planes in (16). Let us consider a 
familiar spacetime slab, call it Q 4 , defined as follows: 

Q 4 = {(t,r)|t a < t < t b , oo < r x , r 2 , r 3 < oo} , 
dQ 4 = {(t,r)|t = t a or t = t b , — oo < r 1 ,r 2 ,r 3 < oo} . (Dl) 
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Next, in classical dynamics we may want to address the following questions: 
Given an ensemble of classical dynamical trajectories {(r(t), p(f)), for t a < t < 
tb}- We take (• • • ) av to represent a classical ensemble average. Then ( i ) what 
is the mean of the difference in r 1 at the final time and initial time? (ii) what 
is the mean of the difference of the squares of r 1 at the initial and final time? 
(Hi) what is the mean of the square of the difference of the final and initial r 1 ? 
and so on. For the ensemble we have the following averaged entities: 


(0 ( rl (^)) av - ( rl ( i «)) av > 

(*0 ( rl (^) 2 ) av - (r\t a ) 2 ) av , 

{Hi) ((r 1 ( t b ) - r 1 (t 0 )) 2 ) av . (D2) 

Let us rewrite (in) as 

(iii) (r'(h) 2 ) m + (r'(C) 2 )„ - (r 1 (tt,)r l (t a )) av - (r'fclr 1 («»))„ 

= (> J (t.)e(f 2 ))„ tt'Slt'.- ( D3 > 

Moving to quantum mechanics, we need to hypothesize, what quantum- 
mechanical entity represents the classical correlation on the rhs of (D3)? If the 
Hamiltonian H is time-independent, a plausible choice is (quantum averages 
with density matrix p given at time t = t a are represented as (• • • )p) 


(m) 


(Ahyfo)). 


1 1 — tb 
ti—t a 


t-2=t-b 

t2—t a 


(D4) 


where the operators are 


r\ta) = r\ 

r x { t b ) = exp(*U(t 6 - t^/tyr 1 exp(-iH{t b - t a )/h). (D5) 


We can also infer that the classical expression 

{iv) ({^{tb) — r ’ 1 (^a)) JV ) av (D6) 

plausibly upgrades into the quantum-mechanical expression 

(») <dt) 

Skipping some steps, we can infer for time-dependent Hamiltonians that (in) 
becomes 




d 3 R 2 


R\R 


1 

2 


x (J 0 l ° 2 (Ti,R i; T2,R 2 ))-] 


T±=tb T2=tb 
Ti—t a T 2 =t a 



(D8) 
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where the last step follows from Gauss’s theorem. 

The lhs of (D8) simplifies when T) = T 2 = t b , say. Since according to (8) and 

(5) the operator U^J°, when integrated over a T=constant slice of spacetime, 
reduces to the unit operator at equal times, this contribution simplifies to a 
single-flux computation: 


J dZR^Rl ) 2 1 )).. 


(38) 


References 

[1] E. Poliak and W. H. Miller, Phys. Rev. Lett. 53, 115 (1984) 

[2] J. Munoz, D. Seidel, and J. G. Muga, Phys. Rev. A 79, 012108 (2009) 

[3] L. I. Schiff, Quantum Mechanics, McGraw-Hill, New York, NY, USA 
(1968), 3rd edition 

[4] L. E. Ballentine, Quantum Mechanics, World Scientific, Singapore (1998). 

[5] G. E. Hahne, Phys. Rev. A 80, 062101 (2009). 

[6] G. E. Hahne, unpublished work. 

[7] G. E. Hahne, J. Phys. A 36, 7149 (2003). An improved version (two 
references were added, but the equation numbers are unchanged) is in 
arXiv:quant-ph/04040 12 . 

[8] D. Stoler and S. Newman, Phys. Lett. A 38, 433 (1972). 

[9] Y. Aharonov and D. Bohm, Phys. Rev. 122, 1649 (1961) 

[10] Y. Aharonov and D. Bohm, Phys. Rev. 134, B1417 (1964) 

[11] I. L. Egusquiza, J. G. Muga, and A. D. Baute, in Time in Quantum 
Mechanics edited by J. G. Muga, R. S. Mayato, and I. L. Egusquiza, 
Springer, Berlin (2008), 2nd edition, Vol. 1, Ch. 10. 


19 



